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While entropy changes are the usual subject of fluctuation theorems, we seek fluctuation relations 
involving time-symmetric quantities, namely observables that do not change sign if the trajectories 
are observed backward in time. We hnd detailed and integral fluctuation relations for the (time 
integrated) difference between entrance rate and escape rate in mesoscopic jump systems. Such 
inflow rate, which is even uirder time reversal, represents the discrete-state equivalent of the phase 
space contraction rate. Indeed, it becomes minus the divergence of forces in the continuum limit to 
overdamped diffusion. This establishes a formal connection between reversible deterministic systems 
and irreversible stochastic ones, conhrming that fluctuation theorems are largely independent of the 
details of the underling dynamics. 

PACS numbers: 05.40.-a, 05.70.Ln 


I. INTRODUCTION 

The past two decades were characterized by the dis¬ 
cussion of fluctuation relations (FRs), or fluctuation the¬ 
orems, valid for systems arbitrarily far from equilibrium 
conditions [ll-[^ (the literature on FRs is vast, see more 
references in the reviews (23 - 1^ ). Initially the focus was 
on deterministic dynamical systems, where Lyapunov ex¬ 
ponents were found to play a fundamental role [ll-U . Al¬ 
though quantities as the phase space contraction rate 
(minus sum of all Lyapunov exponents) are related to 
the entropy production , this concept acquires a more 
immediate physical interpretation in terms of heat flows 
when stochastic systems are considered [§-[l^ [isl - l^ . 
For this reason, and because of the relevance for modern 
technological applications, nowadays it is more frequent 
to discuss FRs for stochastic dynamics. 

The production of entropy, 8 ( 10 ), being related to 
fluxes of heat, matter, etc., changes sign if one goes 
through the trajectory uj backward in time. By now it 
is essentially understood that any FR under examination 
involves a form of entropy production, and that it yields 
a measure of the statistical asymmetry of physical pro¬ 
cesses in time. For example, the integral FR {e~^) = 1 
((...) denotes a statistical average) and the convexity of 
the exponential function provide a statistical mechanical 
derivation of the second law of thermodynamics, (S) > 0. 
Such asymmetry in time exists already at the level of 
trajectories. Entropy production is the physical quantity 
that always determines the time-antisymmetric sector of 
path probabilities [H, [ 2 ^ , 

P{u}) - 

which is generally written as an exponential of some func¬ 
tional A. This includes also a time symmetric (TS) com- 


* baiesi@pd.infn.it 
t falasco@itp.uni-leipzig.de 


ponent K(uj). Yet, such a quantity is completely irrel¬ 
evant in the standard procedure used to prove FRs, in 
which one compares path weights dll) with those of tra¬ 
jectories reversed in time. 

However, it is important to sharpen our understating 
of the meaning and the statistics of TS quantities, since 
it is becoming clear that they are necessary in the for¬ 
mulation of nonequilibrium statistical mechanics. A TS 
quantity that was recognized as a major player in char- 
acterizin g t he evolution of glassy systems is dynamical 
activity j26l - l33l |. which is just the number of jumps (or 
changes of state) that take place during a trajectory, re¬ 
gardless of their direction. Moreover, in the context of 
linear response theory, one finds fluctuation-dissipation 
relations where K{aj) complement the entropy produc¬ 
tion in determining the response of systems far from equi¬ 
librium ( 334 ^ . Mostly these TS observables were con¬ 
tinued to be called dynamical activity, although some¬ 
times also the names “traffic” or “frenesy” |36l | were 
adopted. There are few exainples of fluctuation sym¬ 
metries for TS quantities . Therefore, to better 

understand the physics of nonequilibrium systems, there 
remains the interest of going deeper in this direction and 
find more variants of FRs for TS observables. 

In this paper we briefly introduce some FRs for TS 
quantities. In discrete jump processes we study the in¬ 
flow rate in a state, expressed as a properly defined en¬ 
trance rate minus the standard escape rate. The deriva¬ 
tion of FRs is based on an artificial “auxiliary” dynam¬ 
ics. A similar idea was recently put forward [2l| to define 
the path-space probability of trajectories reversed in time 
when some transitions are only one way (i.e. their reverse 
transition does not exists). With such an approach one 
may draw a generalized FR [2lj. By taking this method 
to the extreme, we will define the auxiliary dynamics as 
that taking place when all jump rates are flipped. It is 
a simple mathematical choice that leads to specific FRs. 
The probability of the auxiliary dynamics, appearing in 
such FRs, can be traced out and one is left with inte- 
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gral FRs, valid for the normal dynamics and in terms of 
sound physical quantities. 

When we perform the limit to overdamped diffusive 
dynamics, we make contact with the approach by Maes 
and van Wieren the auxiliary dynamics here cor¬ 

responds to flipping the force signs and the inflow rate 
becomes minus the divergence of the forces, which is the 
phase space contraction rate of the associated noiseless 
dynamics. In our case the contraction rate is even under 
time reversal, differently from the case of reversible de¬ 
terministic dynamical systems. Nevertheless, an integral 
fluctuation for this quantity is derived. Note that the 
contraction rate can be found, for conservative forces, 
within the definition of the so-called effective potential, 
which is used in computations of reaction pathways [40j| 
and in evaluations of the system activity [4l|, . 


II. JUMP PROCESS 


We consider a system with discrete states {C} whose 
probabilities pt evolve according to the master equation 

dtptiC) = Y, [pt{C')k{C' ^ C) - pt{C)k{C ^ C')], (2) 

C'/C 

where k{C —^ C') is the transition rate from configuration 
C to C. The last term in the master equation contains 
the probability pt(C) times the escape rate from C 

X{C) = Y HCC'). (3) 

C'#C 

Here, next to this concept, we find it useful to define the 
entrance rate e(C) as the sum of transition rates from C 
to C, 


portional to 

n—1 

P{oj) n ^ 

2^0 

n—1 

i=0 


(we are omitting a time-discretization prefactor that is 
common to all trajectories with the same number of 
jumps (3). The standard FR of the entropy production 
is obtained already at the level of single trajectories by 
comparing P(w) with the probability of the path reversed 
in time, which we obtain by applying an involution 0 that 
reverses the order of times (dw = {C{t — s)|0 < s < t}, 
i.e. the initial time of 0 ui is what it was the final time for 
w, etc.). Thus, 


with 


^ pStotiuj) 

P{0UJ) 


( 7 ) 


*5'tot(^) 


n—1 


In pt(C„) + Inpo(Co) -I- P VV(Ci 

i=0 

In pt(C„) -I-In po(Co) + S'(w), (8) 


where S{u}) is the entropy increase in the heat bath. From 
(l7l) one readily observes that an increase of the entropy in 
the bath is associated with processes more likely to take 
place in the normal direction of time, because typically 
> 1 in this case (excluding effects from the bound¬ 
ary terms). Note that by construction this equation picks 
up the time-antisymmetric portion of the path measures 
and completely forgets about the integral of escape rates. 

To obtain a new form of FR we do not simply consider 
time-reversal but we rather define an auxiliary dynamics 
where all rates are replaced by the rates of the inverse 
transitions, 


e{C) = Y HC ^ C). (4) 

C'^C 


In physics often one considers transition rates obeying 
local detailed balance [i^ at temperature T = 1//3, 


fc(C ^ C') 


exp 


fw(c 


C') 


( 5 ) 


where W(C —)■ C) is the work done by the system on 
its environment (a heat bath in equilibrium). The term 
PW{C —>■ C') thus represents the increase of entropy in 
the environment associated to the jump C —)■ C' in the 
system. 

A trajectory from time to = 0 to time t is the time- 
ordered sequence uj = {C(s)|0 < s < t} = {Co = 
C(0),Ci,C2, ■ ■ ■ ,Cn = C{t)} with n jumps Ci_i Ci tak¬ 
ing place at times ti. For a system with initial density 
Pq{C), the probability to observe a trajectory w is pro- 


k*{C ^C') = k{C' ^C), (9) 

so that the “auxiliary” escape rates correspond to the 
entrance rates of the normal dynamics, 

A*(C) = ^ r(C^C') = e(C). (10) 

C'^C 

Note that other choices for “adjoint dinamics” have al¬ 
ready been discussed [IMl- However, those are useful 
to single out specific entropy production terms only. 

The state density p* with the auxiliary dynamics is 
chosen to be the same p we have with the normal dy¬ 
namics. Thus, the path probability under the auxiliary 
dynamics is 

n—1 

p*{uj) ^ pl{Co)e- r(C, ^C,+i) 

2^0 
n— 1 

~ Po(Co)e-7o Y[ fc(C,+i ^ C,), (II) 
2=0 
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while for the corresponding time-reversed path it reads 

n—1 

P*{0u;) ~ l[ k{C^ ^ C.+i). (12) 

i=0 

The initial densities of P*{uj) and P*{9uj) are taken to 
be the initial and final densities of the physical trajec¬ 
tory, respectively, po and pt- Any other choice is equally 
allowed, though (see Sec. IB- 

It is clear that the ratio of the path measure m or 
(fT^ with ® now yields the exponential of a novel path- 
dependent term 

Y{u;)= fdt'[s{Cit'))-\{Cit'))] 

Jo 

d^'7^(C(^')), (13) 

in which there plays a crucial role the instantaneous in¬ 
flow rate, i.e. the imbalance between the entrance rate (|4]) 
and the escape rate dS]) 



n{C) = e{C) - A(C). 
Specifically, the ratio of Eqs. dm) and ( 
= -S{uj)-Yiuj) 

P{UJ) 

while the ratio of Eqs. dni) and is 


IS 


P*{9uj) 

~pW 


— gln[pt(C.i)/po(Co)]-V(a;) 


(14) 

(15) 

(16) 


Using the fact that P{9uj) can be normalized to 1 
(f 'DujP*{9uj) = 1 = (1)*) and the conversion of statisti¬ 
cal averages (...)* = . .P*/P), from (fTHl) we obtain the 

integral FR 


gln[pt(C„)/po(Co)]-V(w)^ _ 


which becomes 

1 for t —>• oo 


(17) 

(18) 


if the boundary contribution is irrelevant. This occurs 
for instance for t —s- oo if the system has a finite number 
of states and is in a stationary regime, pt = po Vt. The 
convexity of the exponential yields also the inequality 
{Y{uj)) > 0. A positive inflow rate is to be expected as 
it is easier on average to jump into the states with high 
p than into those with low p. 

With a similar procedure but using (I15II rather than 
(fT6l) one obtains another FR 


= 1, (19) 


where the entropy production S reappears next to Y. 
Note that (fT7)l and (fT^ a-re valid also in transient con¬ 
ditions. An equation as (fTOl) might be particularly inter¬ 
esting in transient regimes because it does not contain 
explicitly the initial and final density of states. Possi¬ 
bly it might help to study glassy systems, which are the 
paradigm of transient dynamics. 


III. OVERDAMPED DIFFUSION 


In order to identify the transformation of a diffusive dy¬ 
namics equivalent to inverting the rates in a jump pro cess 
we start from the generic Fokker-Planck equation for 
the density pt{x), 

dtPtix) = -dx{pF{x)pt{x)) + Ddlpt{x), (20) 


where pis a mobility, D is a diffusion constant, and F{x) 
is a force. We discretize phase space in tiny units of size 
6, so (EOl) turns into a master equation for configurations 
C = ... — 5,x^x + S,.... With the standard assump¬ 

tion that the dynamics is performed by random walks 
with jumps to nearest neighbors x ^ y = x F S, the 
discretization yields 

dtpt{x) = --^[Fix + 5)pt{x + S) - F{x - S)pt{x - 5)] 
+ ^[pt{x + S) + ptix - S)-2pt{x)], (21) 

that can be rewritten as a master equation with transi¬ 
tion rates 

{ ¥ ~ if 2/ = X - <5 

Hx ^ 2 /) = <^ ( 22 ) 

[ ^ if 2/ = X -b (5. 

When we apply the transformation k*{x y) = k{y —> 
x) we obtain 


D 


fc*(x —>• 2 /) = < 


,52 

D 


S2 


I^F(x) iiy = x-S 
^F(x) iiy = x + 6. 


(23) 


The difference in the state where the forces are evaluated 
is of order 0(1) [1^, and so vanishes when we go back 
to the continuous limit S ^ 0 {D is assumed to be a 
constant). We conclude that changing the sign of F (or 
alternatively of p) gives the transformed rates in terms of 
the original ones, namely kp{x y) = k-p(x —)■ y). The 
path weight associated to dini), describing overdamped 
diffusive dynamics, is 


P{w) ~ po(x(0))exp 


Jo 40 

^^‘dt'5.F(t')| (24) 


and it is straightforward to obtain the equivalent of (El) 
by applying time reversal together with the transforma¬ 
tion F* = —F (or equivalently p* = —p), which again 
gives the probability ratio 


P*{9uj) 


P*tix{t)) _y(^) 

Poixm 


— gln[Pt(a:(t))/po(a:(0))]-V(w)^ 

( 25 ) 
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where now the integral Y is 

Y(uj) = —/i f dt'dxF{t'). (26) 

Jo 

In we found the only previous example where one 
force was flipped to get a FR for TS quantities. Alterna¬ 
tively, we can get again 


Pico) 


(27) 


where the entropy increase in the heat bath is here de¬ 
fined by 


pt 

S{uj) = -^ / dt'F{t')x{t'). (28) 

^ Jo 

Therefore, we find Eqs. (nzi), (HU) and dm) to be valid 
also in overdainped diffusing systems. Besides, from (l2^ 
one can as well derive detailed FRs that link the statistics 
of observables (odd under the joint inversion of time and 
forces) in two systems subject to opposite forces. Exam¬ 
ples of practical interest are systems of non-interacting 
particles exposed to controllable external fields. 

In more than one dimension and with /i = 1, Eq. (USD 
reads 

Y{UJ) = - [ dt'dxMx{t')). (29) 

Jo 

Einstein notation is used from here onward and the vector 
state is denoted as a; = {a;^}. For overdamped stochastic 
systems the inflow rate is thus represented by a function 
with the structure of a divergence of forces, 

7^(a^) = -dxMx) (30) 

when all mobilities are equal to 1. This form has an 
analogous version in deterministic evolution, as discussed 
in the next section. 


IV. ANALOGY WITH THE CONTRACTION 
RATE 


Consider a dynamical systems composed of i = 
1,..., iV degrees of freedom evolving according to deter¬ 
ministic equations 


±i=Fi{x). (31) 

In phase space the system is described by a density p 
which evolves according to the continuity equation, en¬ 
forcing probability conservation: 

dtpix{t),t) = 

=-^xAF^{x{t))p{x{t),t)]. (32) 


Rearranging the terms and introducing the Lagrangian 
derivative d/dt = dt + Xidx^ accounting for time varia¬ 
tions along trajectories, we find 

^\np{x{t),t) = -dxiFAx(t)) =a{x{t)). (33) 

In (1551) we introduced the instantaneous contraction rate 
<7{x), which measures the logarithmic rate of contrac¬ 
tion of phase space volumes —it is identically zero 

in Hamiltonian systems and on average positive in dis¬ 
sipative ones [i^. This is exactly the quantity (1501) for 
which we have derived a FR for the diffusive dynamics. 
Moreover, since the transformation used to find the FR 
in the discrete dynamics (see Sec. El is analogous to that 
used to derive the FR in the diffusive dynamics, we con¬ 
clude that the inflow rate TZ{C) is the equivalent of the 
contraction rate a{x). 

There is indeed a procedure that illustrates this cor¬ 
respondence. In a discrete state system, we spread uni¬ 
formly the occupation probability on the configuration C 
and on its neighbors, i.e. those C such that k{C —>■ C) A 
0. At t = 0 we thus have po{C) = po{C') = 1/V(C), where 
V(C) is the number, or “volume”, of the states centered 
around C. Shortly after the preparation, the time deriva¬ 
tive of the Boltzmann entropy associated to this locally 
flat density is expected to give the logarithmic variation 
of the inverse volume centered on C. Therefore we have 


dt\npt{Cmt=o = V{C)dtPt{C)\t=o 

= V{C)Y,[pt{C')k{C' ^ C) - pt{C)k[C ^ C')] 

c 




= ^ C) - k{C ^ C')] 

c 

= nc), (34) 


which shows that the inflow rate measures the “volume” 
contraction rate. 

The phase-space contraction rate is a key ingredient in 
the derivation of FRs for deterministic dynamics [H, . 

Hence we deem it interesting to strengthen the correspon¬ 
dence between thermostated deterministic systems (i.e., 
Hamiltonian systems with an added nonlinear friction 
term) and stochastic ones applying the ideas of Sec. HIl 
and Sec. m to underdamped diffusion, where the phase 
space does not simply reduce to the configuration space. 
For simplicity we consider two conjugated degrees of free¬ 
dom only, X = {q,mv}, exemplifying a particle of mass 
m and friction coefficient 7 moving in a force field F. 
The motion is described by the Langevin equations 

q = v, mi! = -jv + F{q) + \/2Dy^, (35) 


where ^ is a standard Gaussian white noise with unit 
variance. Since the associated path weight is 


P{uj) ~po(a:(0),u(0))e2* 


>0 




( 36 ) 






5 


the time-integrated contraction rate of the noiseless dy¬ 
namics, Y(uj) = "ft, is singled out by comparison with an 
auxiliary dynamics having negative friction coefficient, 
i.e. 7 * = — 7 . Recalling that time reversal here im¬ 
plies 9v{t') = — the analogous of (IMl) for under¬ 

damped diffusion is 

^ — pln[Pt (»:(*).-«(i))/po(a:(0),«(0))]-7t /'o-r'i 

P(cc) 


of the A^-volume is obtained as 



which gives the time-averaged contraction rate in terms 
of the negative sum of all Lyapunov exponents, 


Essentially, this auxiliary dynamics emulates the time- 
reversal transformation of thermostated Hamiltonian sys¬ 
tems. There [d^, the friction coefficient 7 is indeed re¬ 
placed by the thermostat multiplier, that is an odd func¬ 
tion of the velocity v and thus changes sign upon time 
inversion. Moreover, this auxiliary dynamics is equiva¬ 
lent to that used in Sec. m for underdamped diffusion, 
since "f = Ijfi and we have already noticed that flipping 
the forces is the same as changing the sign of the mobility. 

Similarly to the deterministic case, equa¬ 
tions (EH) and Ell) can be turned into an integral 
FR for the dissipation function H [ 2 ^, 


f dt'rt{x{t')) =\npo{x{0))-\npo{x{t))+ f dt'a{x{t')), 
Jo Jo 

exploiting the freedom in choosing the initial density of 
the backward auxiliary trajectory. Specifically, taking 
= pq in (1251) we obtain the nonequilibrium partition 
identity, i.e. the integral version of the FR 0 



dt'n{t') 


= 1, 


(38) 


valid in overdamped diffusing systems for arbitrary ini¬ 
tial densities and nonconservative forces. The same re¬ 
lation is obtained for underdamped diffusion from (l37l) . 
provided that po is independent of the sign of v, as it is, 
e.g., at equilibrium. 

Finally, it is worth noting that the quantity (1^^ can 
also be related to the (finite-time) Lyapunov exponents 
Ai of the system. The latter are defined by considering 
the growth rates of fc-dimensional volumes supported by 
k linearly independent perturbations Sx^J, with 1 ^ fc ^ 

N [ 13 ]: 


Volfc({5a;W(t)}) ~ . (39) 

The time evolution of such a perturbation 6x in the initial 
condition of (1311) is 

S±i = Fi{x + Sx) - Fi{x) Ki {dx,Fj)5xj = (yF)^jSxj, 

(40) 

which has the solution 

(<) = exp f VF{t')dt'\ 5xj{0). (41) 

\Jo J ij 

Therefore, using the formula relating the determinant to 
the trace, det(expH) = exp(TrH), the relative variation 



N 


y:A,. 


(43) 


If (1^ is stochastic, as in the overdamped diffusion 
considered in the previous section, the Lyapunov ex¬ 
ponents identify the separation rate of nearby trajec¬ 
tories subjected to the same noise realization, and 
Eqs. (1251) and EH are statements about the statistics 
of their spectrum. Note that Eqs. EH) and ( 1 ^ can be 
restricted to quantities measured on subspaces of the full 
phase space using auxiliary dynamics where only some 
components of the forces are changed in sign. In partic¬ 
ular, when F* = —Fi if i = k and F* = Fi otherwise, we 
get for example (no Einstein summation) 



1 

15 


dt' Fk{t')xk(t') 


dt'd^,Fk{t') 


= 1 . 
(44) 


Since a fc-volume is obtained as the determinant of the 
k X k minor of 5xi{t)/5xj[d), Eq. (l44l) relates a single 
Lyapunov exponent Aj,, 


= Volfc(t) - Volfe_i(t) = exp dt'dx^Fkij')^ , 

(45) 

to the entropy produced in the k direction. 


V. CONCLUSIONS 

We have introduced the concept of inflow rate. The 
fluctuations of its time integral are characterized by two 
fluctuation relations, one involving also boundary terms 
(density of states) and one involving the entropy produc¬ 
tion. In the continuum limit bringing jump systems to 
overdamped diffusion, the inflow rate becomes (minus) 
the divergence of forces, i.e. the phase space contraction 
rate. This quantity in reversible deterministic dynamical 
systems is known to be odd under time reversal. How¬ 
ever, in general the phase space contraction rate can have 
any or no parity. In particular, in our case it results 
even under time reversal, due to the fact that the over¬ 
damped Langevin dynamics is not time reversible. Note 
that time reversibility is a condition commonly invoked 
in the derivations of FRs, hence it is sufficient, but it is 
known to be unnecessary [sil - f^ . 

We hope that these findings will be useful in devel¬ 
oping nonequilibrium statistical mechanics —e.g. in de¬ 
riving symmetry properties of response coefficients far 
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from equilibrium, as much as standard FRs are employed even parity under time reversal, 
around equilibrium to obtain the Onsager reciprocity re¬ 
lations |54| — for which by now it is clear that we need 

not only to have entropy flows under control, but also to ACKNOWLEDGMENTS 

better understand the statistics of activities that have an 
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